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Forgetting in the Synchronization of Quantum Networks 


Shuangshuang Fu, Guodong Shi, and Ian R. Petersen* 


Abstract 

In this paper, we study the decoherence property of synchronization master equation for networks of 
qubits interconnected by swapping operators. The network Hamiltonian is assumed to be diagonal with 
different entries so that it might not be commutative with the swapping operators. We prove a theorem 
establishing a general condition under which almost complete decohernece is achieved, i.e., all but two 
of the off-diagonal entries of the network density operator asymptotically tend to zero. This result 
explicitly shows that quantum dissipation networks tend to forget the information initially encoded 
when the internal (induced by network Hamiltonian) and external (induced by swapping operators) 
qubit interactions do not comply with each other. 

Keywords: quantum networks, synchronization, decoherence 

1 Introduction 

Inspired by the developments of distributed consensus control for classical network systems in the past 
decade [DEI El, consensus and synchronization problems of quantum networks have also recently attracted 
attention in the research community mmnm- Sepulchre et al. [6] generalized consensus algorithms to 
non-commutative spaces and presented convergence results for quantum stochastic maps, and showed 
how the Birkhoff theorem can be used to analyze the asymptotic convergence of a quantum system to a 
fully mixed state. Mazzarella et al. [7] made a systematic study regarding consensus-seeking in quantum 
networks, introducing several classes of consensus quantum states and a quantum generalization to the 
gossip iteration algorithm based on pairwise swapping operators for reaching a symmetric state (consensus) 
over a quantum network. The class of quantum gossip algorithms was further extended to symmetrization 
problems in a group-theoretic framework [8j. 

The analysis of quantum consensus seeking was further developed using the graphical methods for 
studying classical network systems [1], and it was shown that the vectorized density operator evolving 
along quantum consensus dynamics is equivalent to a number of parallel classical consensus dynamics 
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over disjoint subgraphs [9], which enabled us to study quantum consensus dynamics via their classical 
analogous with all details inherited. Furthermore, when the network Hamiltonian is commutative with 
the swapping operators, one can derive a so-called quantum synchronization master equation [9] as the 
quantum counter part of the classical linear synchronization results H3H5I . This quantum synchronization 
master equation can be physically realized via quantum dissipation networks where quantum nodes are 
interconnected by local environments un¬ 
it however has been understood that when a quantum system interacts with the environment through 
dissipative couplings, the quantum information encoded in the system is often washed out in the sense that 
the off-diagonal entries of the system density operator asymptotically vanish. This phenomenon is known 
as decoherence. In this paper, we study the decoherence property of synchronization master equation for 
quantum networks of qubits interconnected by swapping operators. The network Hamiltonian is assumed 
to be diagonal but with different diagonal entries so it might not be commutative with the swapping 
operators. We prove a theorem establishing a general condition under which almost complete decohernece 
is achieved, i.e., all but two of the off-diagonal entries of the network density operator asymptotically tend 
to zero. This result explicitly shows that quantum dissipation networks tend to forget the information 
initially encoded when the internal and external qubit interactions, respectively induced by the network 
Hamiltonian and the swapping operators, do not comply with each other. 

The remainder of the paper is organized as follows. Section [2] introduces the problem definition and 
presents the main result obtained. A brief introduction to the quantum mechanics related to the devel¬ 
opments of the current paper is also provided in Section [2] as well as a few numerical verifications of the 
theoretical result. Section [3] establishes the detailed proof of the main result, where the arguments are 
organized into step-by-step blocks. Finally Section [4] concludes the paper. 

2 Problem Definition, Main Result, and Examples 

In this section, we define the problem of interest, present the main result, and provide numerical examples 
illustrating the obtained result. 

2.1 Quantum States, Density Operators, and Partial Trace 

We first give a brief introduction to quantum systems’ states. We refer the readers to [5] for a comprehensive 
treatment. 

2.1.1 Quantum States 

The state space associated with any isolated quantum system is a complex vector space with inner product, 
i.e., a Hilbert space H. The system is completely described by its state vector, which is a unit vector in 
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the system’s state space and often denoted by \ip) G H (known as the Dirac notion). The state space of a 
composite quantum system is the tensor product of the state space of each component system, e.g., two 
quantum systems with state spaces Ha and Hb, respectively, form a composite system with state space 
Ha ®Hb, where <g> stands for tensor product. If the two quantum systems are isolated respectively with 
states IV’a) G Ha and \^b) G Hb, the composite system admits a state \if a) <8> |Vib)- 

2.1.2 Density Operators 

For an open quantum system, its state can also be described by a positive (i.e., positive semi-definite) 
Hermitian density operator p satisfying tr(p) = 1. A quantum state \ip) E H , induces a linear operator, 
denoted by 

lY^lfV)) = (|V’),|*))| 

with being the inner product^ equipped by the Hilbert space H. Then p = |'i/’)(V ; l defines the 

corresponding density operator. Density operators provide a convenient description of ensembles of pure 
state : If a quantum system is in state | i/>i) with probability pi where )T) ?: Pi = 1) its density operator is 

p = 

i 

Any positive and Hermitian operator with trace one defines a proper density operator describing certain 
quantum state, and vice versa. 

2.1.3 Partial Trace 

Let Ha and Hb be the state spaces of two quantum systems A and B, respectively. Their composite 
system is described by a density operator p AB . Let 2a, 2b, and 2ab be the spaces of (linear) operators 
over Ha, Hb, and Ha ®Hb, respectively. Then the partial trace over system B, denoted by Tr% s , is an 
operator mapping 2ab to 2a defined by 

^H B (\PA){qA\ <S> |pb>(<?b|) = |PA)(9A|Tr(jp B )(g B |) 

for all | pa), | <?a) G Ha, | Pb), \qb ) G H B - 

The reduced density operator (state) for system A , when the composite system is in the state p AB , is 
defined as p A = Tr% B (p AB ). The physical interpretation of p A is that p A holds the full information of 

. A D 

system A m p . 

1 Under Dirac notion this inner product is written as (| if>), |a;)) = (il>\x), where (ip\ is the dual vector of | ip). 
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2.2 Qubit Network and Its Synchronization 

In quantum systems, the two-dimensional Hilbert space forms the state-space of qubits (the most basic 
quantum system). Let LI be a two-dimensional Hilbert space for qubits. The standard computational basis 
of is denote by |0) and |1). An n-qubits quantum network is the composite quantum system of n qubits 
in the set V = {1,..., n}, whose state space is the Hilbert space 'H <Zn = T~t 0 ■ ■ ■ 0 T~L, where 0 denotes the 
tensor product. The swapping operator between qubits i and j, denoted as Uij , is defined by 

Uij(\qi) 0 • • ■ 0 \qi) 0 ■ ■ ■ 0 \qj) 0 ■ ■ ■ 0 |g n )) = 
ki) 0 • • • 0 \qj) 0 • • • 0 | qi) 0 ■ ■ ■ 0 | q n ), 


for all qi E {0,1 },i = 1 ,,n. In other words, the swapping operator Uij switches the information held 
in qubits i and j without changing the states of other qubits. 

The density operator of the n-qubit network is denoted as p. A quantum interaction graph over the 
n-qubit network is an undirected, connected graph G = (V, E), where each element in E, called a quantum 
edge, is an unordered pair of two distinct qubits denoted as {i, j} E E with i,j E V. The state evolution 
of the quantum network is given by the following master equation [9], 

t = ~i lH ’ p]+ £ ( u n<p u k-p)- f 1 ) 

LAKE 

where [•, ■] denotes the commutator of two operators, H is the effective Hamiltonian as a Hermitian operator 
over the underlying Hilbert space, i 2 = — 1, h is the reduced Planck constant, Ujk is the swapping operator 
between j and k. As discussed in [9], the above synchronization dynamics is a Markovian master equation 
in the Lindblad form OEM and can be physically realized via building suitable local environments among 
the qubits [Ill- 

Let P be the n’th permutation group and assume the initial time is 0 for the system (jTj) . It has been 
shown in [9j that when the network Hamiltonian H is commutative with the swapping operators, i.e., 
[H, Ujk] = 0 for all {j, k} E E, quantum synchronization is achieved in the sense that (cf., [9]) 


hm (p(i) - e~ lHt / h p*e lHt / h '\ = 0 (2) 

along the system 

qubit k at time t. For the limiting trajectory, there holds for all j = 1 ,,n that 


, where p* = ^t Zk e p U^p{Qi)Ul. Let p k (t) := Trbe the reduced state of 


TV, 




which in turn leads to 


(3) 


^lim (p k {t) — p m (t)^j =0, k, m E V. 


(4) 
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2.3 Main Result: A Quantum Forgetting Theorem 

Throughout the remainder of the paper, we investigate p(t) under the following standard basis of 7i &n : 

B := j|<?i... q n )('Pi ■ ■ ■ p n I : Pi,Qi G {0, l},i G N/j. 

We identify the operators with their matrix representations under the basis B, for the ease of presentation. 
We denote [p(t)] ^ ? ^ p | as the \q\ ■ ■ ■ q n )(pi ■ • -Pn|-entry of the density operator p(t) under the basis 
B. The diagonal entries of the elements in B are put in the set 

®D := jbi • • -Pn){Pl • ■ ■ Pn | • Pi G {0,1 },i G V}. 

For simplicity we always write z = z±... z n with Zi G {0,1}, i G V. In this paper, we are interested in the 
evolution of the system 0 in the absence of the commuting condition between the Hamiltonian and the 
swapping operators. Particularly, we are interested in the decoherence of the system i-e., decaying of 
the off-diagonal entries of the density operators. To be precise, we introduce the definition of decoherence 
in the following. 

Definition 1 The system achieves \x)(y\-decoherence for if [p{t)] 1 = 0. 

We impose a standing assumption on the network Hamiltonian H. 

Assumption There are 2 n real numbers A| p ^ p | G R, \p)(p\ G Bd such that H = X^|p)(p|eB D A| p )< p | \p)(p\ ■ 

Under the above assumption, the network Hamiltonian H is diagonal under the standard basis. Since the 
Hamiltonian H is a Hermitian operator, one can always find a basis of Tl® n under which H is represented 
by a diagonal matrix. The assumption that H is diagonal under the standard basis is however quite 
restrictive. Nevertheless, this assumption allows for basic non-commuting properties between H and the 
swapping operators, and in the meantime enables us to derive some explicit result for the decoherence of 
the system 0. 

Let C,k be the combinatorial number of selecting k from n objectives. Denote 0 = 0 ... 0 and 1 = 1... 1 
both with n digits. The following is our main result. 

Theorem 1 The following statements hold for the system 0 - 

(i) If the elements A| p ^ p | — A| p /^ p /|, p p' are pairwise distinct, then almost complete decoherence is 
achieved in the sense that \x)(y\-decoherence is reached for all x ^ y satisfying either x ^ {0,1} or 
y {0,1}, for which the convergence is at an exponential rate; 

(n) [p(*)]|o><i| = [p( 0)]|o) ( i|e _ * (A l°><°l _A l 1 >[p(f)]| 1>(0 | = [p( 0 )]| 1)(0 |e-< A ii>(irAo)(o|Wn. 
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(Hi) For x = x\ ... x n with Xi = k, there holds that 


A™ W‘)]|*)<*i = E W°)]|») ( »|/ C " 

y : J27= i Vi= k 

where the convergence is also exponential. Consequently, there are at most n + 1 different values for 
the limits of the diagonal entries of p(t ). 


Note that there are only two off-diagonal entries, [p(t)]|o)(i| and [p(t)]|i)(o|> that can possibly be not 
vanishing for the system |lj). It is clear that if the elements A| p )( p | — A| p q( p /|, p / p' are pairwise distinct, 
the network Hamiltonian is no longer commutative with the swapping operators in the system 0 The 
decoherence result established Theorem [I] reveals that the quantum network along the system ([Tj) then 
tends to forget almost all the information contained in the off-diagonal entries of the initial value p( 0 ), 
which represent initial correlations among the the basis states [3]. Furthermore, if the condition that the 
elements A| p \/ p | — A| p /w p /|, p / p’ are pairwise distinct does not hold strictly, it is clear from the proof 
of Theorem [l] that a network Hamiltonian H being non-commutative with the swapping operators in the 
system 0 continues to tend to wash out the off-diagonal entries of the network density operator, just 
possibly leaving a few nonzero off-diagonal entries. 


2.4 Numerical Example 


In this subsection, we present a numerical example to illustrate the obtained main result. We consider 
three qubits indexed in V = {1,2,3}. Their interaction graph is fixed as the complete graph, i.e., E = 
{{1, 2}, {2,3}, {1, 3}}. Let a \2 = «i 3 = «23 = 1- We denote x = X\X 2 X 3 with Xi £ {0,1}, and whenever 
applicable we identify x as a binary number. The initial network state is chosen to be 


Po = (X]l x ))(lZ^l)/ 128+ ( ^ (x + l)|s>(x|)/72. 

* x |:r)( 3 :|eBD 


The network Hamiltonian is chosen to be 


H= J2 ^ X \x){x\ 

|s)(®|eB D 

so that our standing assumption is satisfied. 

We first plot the evolution of [p(0] | x \( x | f° r x . Clearly the eight trajectories are asymptotically 
grouped into four clusters. We also introduce 


Eo(t) := E 

\x)(y\: ^y,|z><^{|0>(l|.|l><0|} 


\x)(y\ 


(5) 


as a measure of decoherence for all off-diagonal entries of p(t) except for |0)(1| and |1)(0|. We also plot 
E 0 (t) and clearly it tends to zero exponentially. 
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Figure 1: The evolution of the diagonal entries of the network density operator (left) and the E a (t ) (right). 

3 Proof of the Main Result 

In this section, we present the proof of Theorem [l] The analysis is based on splitting the entries of the 
density operator into decoupled subgroups where interactions only take place inside each subgroup. The 
idea of breaking down large density operators of multiple qubits can in fact be traced back to [10] using 
Stokes tensors. In [9], the method of investigating the interconnection of the individual entries of the 
network density operator was systematically studied. 

3.1 Graphical Decomposition 

We first establish a technical lemma. 


Lemma 1 Denote C as an operator over the space of p by C (p) := [H,p\. Then [Cjuw^i = (Auw-eI — 

X \y)(y\)\ x )(y\- 

Proof. Based on our standing assumption on the definition of H, we obtain 

[^-'1 \x) {y\ = ^ X \p)(p\\p)(p\ x ><y|- E \p) (p\\ x ) (y\p) {p\ 

|p>(p|sb d b>(p|eB D 

= X] x \ P )(p\ 5 (p,x)\p)(y\ - x \p){ P \ 6 (Piy)\ x )(p\ 

|p>(p|sb d Ip>(p|gb d 

— (A|x) (:e| ^ll/Xl/l) I 3 ') (y\ 1 (®) 

where 6(a, b) = 1 if a = b and 5(a, b) = 0 otherwise. This completes the proof. □ 

We also recall the following lemma, which is a variation of the Lemma 4 in [9J. 

Lemma 2 Let Uij in the swapping between i and j in the permutation group P. i.e., Uij(i ) = j, Uij(j) = i, 
and Uij{k ) = k for k ^ i,j. Then there holds Ujk\x)(y\Uj k = \ujk(x))(ujk{y)\- 
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In light of Lemmas [I] and [2j we can now rewrite the system 0 into its entry-wise equivalence: 

dt \x){y\ — ^IwXs/l) [x(^)] \x){y\ + ([xW] \u jk (x)){u jk (y)\ [p^)] | x )({/|) ’ ^ 

{j,k} SE 

where x = x\ ... x n , y = y\ . . . y n with Xj, y* G {0,1}, i G V. We see from Q that (cf., Lemma 5 [9]) 

^ \xi...Xn)(y\...y n \ := ■ ■ • x ir(n)) ■ ■ ■ 2/7r(n) I • ^ ^ 

forms a subset of entries whose state evolution is not influenced by entries outside. It is clear that if either 
x (f {0,1} or y (jz {0,1} holds, then |<y|I > 2 from Lemma [5] of Moreover, Theorem JTJ (ii) follows 
from direct calculation from the system (|7j) since 7H|o)(i| an d 7^|i)(o| are singletons. 

We are now ready to state the following lemma which transforms the decoherence of system (jT]) to a 
synchronization problem of a classical network. 

Lemma 3 Consider a classical network with N nodes indexed in the set V = {1,... ,N} with an under¬ 
lying interaction graph G = (V, E) which is undirected and connected. Let node i possess a state G C. 
The evolution of the Xi is given by 

^Xi(t) = iBMt) + J2 fa® ~ Xt(tj) (8) 

j:{i,j}eE 

where 6i G IR for all i G V. The following statements are equivalent. 

(i) The system 0 satisfies lirn^oo Xjft) = 0 for alii G V if N > 2 and the 8i, i G V are pairwise distinct; 

(ii) The system |7p achieves \x)(y\-decoherence if \R\ x )(y \| > 2 and the A| p ^ p | — A| p q( p /|, p / p' are pairwise 
distinct. 

Proof. Denoting 6; = — (A^^i — \ y ){ y \) /h and investigating the system 0 over the set 1Z\ x )(y\ with x / y, 
the desired equivalence becomes clear from the definition of decoherence immediately. □ 

3.2 A Classical Detour 

We proceed to make a further investigation to the system ([8]). To this end, we make use of the realification 
method to investigate the system ([8]) via studying the real and imaginary parts separately. We write 

Xi(t) = Ri(t) + iSi(t) 

where Ri(t) and Si(t) are the real and imaginary components of Xf(t), respectively. Denote Yi(t) = 
( Ri(t ) Si(t)) T for i G V. Then the system Q reads 

jY^t) = AiYi{t) + (Xj(t) — Yi(t)), i G V 
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(9) 


where 


Ai = 




i 6 V. 


Clearly the system defines a classical linear synchronization problem with non-identical node self 
dynamics specified by the Ai (cf., [Q, T5j). The following is an intermediate result for the system (|8]) 
established by studying its realification system Q. 


Lemma 4 Denote f(t ) := maxjgy ||X;(t)|| 2 . Then f(t) is a non-increasing function along the system §. 


Proof. Clearly /(f) is a continuous but not necessarily continuously differentiable function. In this step, 
we prove that /(f) is a non-increasing function along the system ®> by showing that its Dini derivative is 
always non-positive. 

The upper Dini derivative of a function h : (a, b) —> 1R at f € (a, b) is defined as [16] 

PV»W = lims U p ' l(t + S) — Mt> . 


s—»0+ 


ii 112 

Define Z(f) := arg maxj g y l](t) . The Lemma 2.2 of [T7] enables us to derive 

D + fit) = max — ||y](t)|| 2 
J K ’ iex(t) clt 11 v ;|1 


2 max (Y i (t),A i Y i (t)+ V (Yj — Yf)\ 
iei(t) \ . —■* / 




a) 


b) 


^ ,T.,0 K ' 


< - max (||Ti(t)|| 2 - ||lj(t)|| 2 ) 


iex(t) 




c) 

< 0 , 


( 10 ) 


( 11 ) 


where a) is based on the fact that (Yi(t), AiYift)) = 0 from the definition of Ai, b ) follows from the 
elementary inequality a T b < (||a|| 2 + ||6|| 2 )/2 for two vectors a and b, and c) is due to the definition of 
Z(t). Based on the properties of the Dini derivative, © leads to that /(f) is a non-increasing function 
along the system ([9]) for all f > 0. □ 

We are now ready to prove the following key lemma for the system ([8]). 


Lemma 5 For the system |<§j) with initial value X (0) = (Xi(0)... Jy(t)) T , there exists a non-negative 
real number Z X ( o) > 0 such that lim t _>.oo ||-Xj(t)|| = Z x to) f or a M * £ V. 

Proof. The analysis will be carried out for the system ([9]). Since /(f) is a non-increasing by Lemma[4j for the 
initial value X (0) = (X\ (0)... X x (t)) T , there exists a constant f*(X( 0)) > 0 such that lim^oo /(f) = /*. 
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We prove the desired lemma by showing that lim^oo ||l^(t)||“ = /* for all i e V via a contradiction 
argument. 

Suppose there is a node *o E V satisfying g* := liminf^oo ||Y^ 0 (i)|| 2 < /*. Consequently, there exists 
an infinite time sequence t\ < ■ ■ ■ < t m < ■ ■ ■ such that 


Yi 0 (t m )\\ 2 < + /*), m = 1,2,.... 


( 12 ) 


On the other hand, from the definition of /* and the analysis of Step 1 we conclude that for any e > 0, 
there exists T e > 0 such that 

\\Y(t)\\ 2 < f* + e, t>T e . (13) 

We build the remainder of the proof in steps. 

Step 1. Take a time instant t m and without loss of generality let t m > T e . In this step, we bound ||Tj 0 (t)|| 2 
during the time interval [t m ,t m + !]• Similar to the derivation of (jTTj) , we have 

d '\Y i0 (t)\\ 2 =2(Y i0 (t), £ {Yj(t)-Y i0 (t))) 

j-{io j}eE 


dt 


< £ (lln(*)|| 2 -||r i0 (t)|| 2 ) 

j-{i odleE 

<{n- 1 )(/* + e — ||y io (f)|| 2 ^ 
for all t >T e . Invoking the Gronwall’s inequality, we further conclude 

ll^oWf < e -("- 1 X*-^)||y io(im) || 2 + (l - (/* + e), t > t r 


Plugging in (12), (15) leads to 


(14) 


(15) 


(16) 


\\Yio{t)\\ < (g* + (l — c) (/* + e ) 5 f 6 [tm,tm + 1]) 

where ( = e^^ n ^^/2. 

Step 2. Now that the graph G is connected, there must be a node i\ / iq such that {io,*i} E E. In this 
step, we bound ||yq(f m + 1)|| 2 . For ||Tj 1 (t)||", we have 
d 


dt j^ 1 w|| 2 <(||^ 0 w || 2 - ||^i w|| 2 ) + E (ll^wil 2 - Uncoil 2 ) 

j^=i odhJjeE 

<(g* + (l — C) (/* + e ) — ||^ii (t) || + (n — 2) ^/* + e — ||li 1 (t) || ^ 


(17) 


for all t E [tm ) tm + 1], where in the second inequality we have used (18). Again, invoking the GronwalPs 
inequality, we conclude 

11^*1 (^m + 1)|| < e ( n 1 '*||^ii(^m)|| + (l — e ^) (^Cd* + (l — C) (/* + e ) + ( n ~ 2)(/* + e)^ /{jl — 1) 

< e ( n 1 )(/* + e) + (l — e 1 ' ) ) ((g* + (l — C) (/* + e ) + ( n ~ 2)(/* + e)')/(n — 1) 


< (<K)3* + (f - (0C)) (/* + e), 


(18) 
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where (f> = (1 — e ( n 1 ^)/(n — 1). In fact, we even know 

|| Y s (t m +1)|| 2 < ( 4>C)g* + (f - (<K))(/* + e), s € {i 0 ,n} (19) 

since tf> E (0,1). 

Step 3. Since the graph G is connected, we can recursively apply the arguments in the Steps 2 and 3 to 
the rest of the nodes, and eventually establish 

+ ti — 1)|| < (<t>C) n l Q* + (l — {4 > C) n 1 )(/* + e )> sEV. (20) 

This implies 


f (t m + n — 1) < (4>C) n l g* + ^1 — (4>() n 1 ^(/* + e), 


( 21 ) 


which contradicts the definition of /* if 


OC) n 

1 - (<(>C) n_1 

Therefore, we have proved that liminf^oo ||Y)(t)|| 2 = /* for all i 
holds limsup^oo ||Y)(t)|| 2 < /* in light of Lemma |ij Consequently, 
/* for all i E V, which completes the proof. 


( 22 ) 

E V. On the other hand, there always 
we have shown that Hindoo ||lq(t)||~ = 

□ 


3.3 Completion of the Proof 

In this subsection, we complete the proof of Theorem [l] 


3.3.1 Decoherence 


With Lemma [3j the decoherence statement for the system ([Tj) holds if lim^oo Xi(t) = 0 for all i E V for 
the system ([8]) with N > 2 when the i E V are pairwise distinct. In fact, we are going to show a slightly 
stronger result for the system ([8]) which only requires that there exist two distinct values within the Oj. 

We recall a few preliminary results on the limiting set of autonomous systems. Consider the following 
autonomous system 

x = f(x), (23) 


where / : R“ —> R“ is a continuous function. Let x(t) be a solution of (23) with initial condition x(to) = ar 


Then Pq C R rf is called a positively invariant set of (23) if, for any to £ R and any a; 0 E Pq, we have 


x(t) E flo) t > to, along every solution x{t) of (23). 

We call y a w-limit point of x(t) if there exists a sequence {t^} with lim^oo tj~ = 00 such that 
lirn^oo x(tfc) = y. The set of all w-limit points of x(t) is called the w-limit set of x(t), and is denoted as 
A+( x(t)j. The following conclusion is well-known na¬ 


il 






Lemma 6 Let x(t) be a solution of ^25 | ). Then A + (x(f)) is positively invariant. Moreover, if x(t) is 
contained in a compact set, then A + (x(t)) ^ 0. 

We are now ready to state the following result for the system Q. 

Proposition 1 For the system 0 with N >2, the following statements hold. 

(i) limt-^oo Xi(t) = 0 for all i 6 V if there exist at least two distinct values within the 6i. 

(ii) lim^oo ||Xj(t) — e ^*Si=i^iM|| = 0 if there is 9 £ R such that 9 = 9i for all i € V. 

Proof, (i) By Lemma[4j for any given initial value, the trajectory of the realification system Q is contained 
in a compact set. Therefore, the uj -limit set of Y(t) = (Y\ (t)... YN{t)) T along the system ([ 9 ]) is nonempty 
and invariant in light of Lemma [BJ On the other hand, making use of the Lemma [5j one finds that the 
w-limit set of Y(t), denoted A + (T(f)), can only be a subset of the set 

A := {Y* = (Yf . ..Yfj) : ||F/|| = Z v(0) }. (24) 

Without loss of generality we assume Z X ( 0 ) > 0 since otherwise the desired result holds immediately. The 
remaining argument relies on showing that any subset of A cannot be invariant for the system ([ 9 ]) if there 
exist at least two distinct values within the 0*. We only need to establish two facts. 

FI) For any T* E A + (T(f)), there must hold Yf = ■ ■ ■ = Yjy. This is due to that as long as Y* ^ Y*, 
the trajectory starting from W must leave the set A since the terms A{Y* are always perpendicular 
to the tangential directions of the manifold A. See Figure [2] for an illustration. 

F2) From FI, we have A+ (Y(t)) C A fj T with T := jV* = (F* ... Y*) : Yf = ■ ■ ■ = Y* }. However, T 
cannot be invariant if there are at least two distinct values within the Oi since Zx( 0 ) > 0. 

Therefore, one must have Z\(o) = 0 anc ^ desired conclusion follows. 

(ii) The conclusion is straightforward using the transformation Xj(t) = e~ lBt Xi(t), where clearly 

jXi{t)= Y, - Xi(t)) 

j:{i,j}£E 

and thus there holds that lim^oo ||Aj(t) — 1| = 0 . 

This completes the proof. □ 

Combining Lemma [ 3 ] and Proposition [IJ Theorem [I] (i) is immediately proved, where the exponential 
rate of convergence is simply resulted from the linear structure of the system 0. 
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Figure 2: Illustration to why there must hold Y* = • • ■ = Yy for any Y* E A + (Y(f)): (i) If Y/* / Y)* and 
there is a link between i and j, the trajectory starting from Y* must leave A since the terms AjY* are 
always perpendicular to the tangential directions of the manifold A; (ii) such a pair always exists if there 
are i,j E V with Y* ^ Y* since the graph G is connected. 

3.3.2 Diagonal Entries 

From the system (J7J), we have 

^[P^\\x)(x\~ (l/^)] \u jk (x)){uj k (x)\ [P^)] |a;> <a;|) ’ 

{j,k} SE 

which is consistent with the case when the network Hamiltonian H is commutative with the swapping 
operators. Theorem [lj (iii) readily follows from the analysis established in |9] by applying Theorem 1 and 
Lemma 5 of [5j. 

4 Conclusions 

We have made a further investigation to the decoherence property of synchronization master equation for 
quantum networks of qubits interconnected by swapping operators. The network Hamiltonian is assumed 
to be diagonal but with different diagonal entries so it might not be commutative with the swapping op¬ 
erators. We proved a theorem establishing a general condition under which almost complete decohernece 
is achieved, i.e., all but two of the off-diagonal entries of the network density operator asymptotically 
tend to zero. This result explicitly revealed that quantum dissipation networks would forget the informa¬ 
tion initially encoded when the internal (network Hamiltonian) and external (swapping operators) qubit 
interactions do not comply with each other. In future, it is interesting to look at the case with switch- 
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ing interactions where nontrivial coherence could be left if the switching signal properly responds to the 
network Hamiltonian. 


References 

[1] J. Tsitsiklis, D. Bertsekas, and M. Athans, “Distributed asynchronous deterministic and stochastic 
gradient optimization algorithms,” IEEE Trans. Autom. Control , vol. 31, no. 9, pp. 803-812, 1986. 

[2] A. Jadbabaie, J. Lin, and A. S. Morse, “Coordination of groups of mobile autonomous agents using 
nearest neighbor rules,” IEEE Trans. Autom. Control , vol. 48, no. 6, pp. 988-1001, 2003. 

[3] R. Olfati-Saber and R. M. Murray, “Consensus problems in networks of agents with switching topology 
and time-delays,” IEEE Trans. Autom. Control , vol. 49, pp. 1520-1533, 2004. 

[4] M. Mesbahi and M. Egerstedt. Graph Theoretic Methods in Multiagent Networks. Princeton University 
Press. 2010. 

[5] M. A. Nielsen, and I. L. Chuang. Quantum Computation and Quantum Information. 10th Edition. 
Cambridge University Press, 2010. 

[6] R. Sepulchre, A. Sarlette and P. Rouchon, “Consensus in non-commutative spaces,” Proc. f9th IEEE 
Conference on Decision and Control , pp. 6596-6601, Atlanta, USA, 15-17 Dec., 2010. 

[7] L. Mazzarella, A. Sarlette, and F. Ticozzi, “Consensus for quantum networks: from symmetry to 
gossip iterations,” IEEE Trans. Autom. Control , 60(1): 158-172, 2015. 

[8] L. Mazzarella, F. Ticozzi and A. Sarlette, “From consensus to robust randomized algorithms: A 
symmetrization approach,” quant-ph, arXiv 1311.3364, 2013. 

[9] G. Shi, D. Dong, I. R. Petersen, and K. H. Johansson, “Reaching a quantum consensus: master equa¬ 
tions that generate symmetrization and synchronization,” IEEE Transactions on Automatic Control , 
in press, arXiv:1403.6387, 2015. 

[10] C. Altafini, “Representing multiqubit unitary evolutions via Stokes tensors,” Physical Review A, vol. 
70, 032331, 2004. 

[11] F. Verstraete, M. M. Wolf, and J. I. Cirac, “Quantum computation and quantum-state engineering 
driven by dissipation,” Nature Physics , vol. 5, pp.633-636, 2009. 

[12] G. Lindblad, “On the generators of quantum dynamical semigroups,” Comm. Math. Phys., vol. 48, 
no. 2, pp. 119-130, 1976. 


14 


[13] H.-P. Breuer and F. Petruccione. The Theory of Open Quantum Systems. Oxford University Press, 
2002 , 1st edn. 

[14] C. W. Wu and L. O. Chua, “Synchronization in an array of linearly coupled dynamical systems,” 
IEEE Trans. Circuits Syst ., vol. 42, pp. 430-447, 1995. 

[15] C. W. Wu. Synchronization in complex networks of nonlinear dynamical systems. World Scientific, 
2007. 

[16] J. Danskin, “The theory of max-min, with applications,” SIAM J. Appl. Math., vol. 14, pp. 641-664, 
1966. 

[17] Z. Lin, B. Francis, and M. Maggiore, “State agreement for continuous-time coupled nonlinear sys¬ 
tems,” SIAM J. Control Optim., vol. 46, pp. 288-307, 2007. 

[18] N. Rouche, P. Habets, and M. Laloy, Stability Theory by Liapunov’s Direct Method. New York: 
Springer-Verlag, 1977. 


15 



